Mixing and segregation rates in sheared granular materials 
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The size-segregation of granular materials, a process colloquially known as the Brazil Nut Effect, 
has generally been thought to proceed faster the greater the size difference of the particles. We 
experimentally investigate sheared bidisperse granular materials as a function of the size ratio of 
the two species, and find that the mixing rate at low confining pressure behaves as expected from 
percolation-based arguments. However, we also observe an anomalous effect for the re-segregation 
rates, wherein particles of both dissimilar and similar sizes segregate more slowly than intermediate 
particle size ratios. Combined with the fact that increasing the confining pressure significantly 
suppresses both mixing and segregation rates of particles of dissimilar size, we propose that the 
anomalous behavior may be attributed to a species-dependent distribution of forces within the 
system. 

PACS numbers: 45.70.Mg, 47.57.Gc, 81.05.Rm, 64.60.ah 



Accurate knowledge of the rate at which granular ma- 
terials segregate by size under shear |TJ [2] is significant 
for such applications as avalanche hazard prediction and 
the design of industrial particle separation chutes. One of 
the most common segregation phenomena, the Brazil Nut 
Effect, is broadly observed and has been associated with 
a variety of proposed mechanisms [3]. For shear flows, ki- 
netic sieving theory in various forms [HOE] has been the 
most promising. These theories rely on statistical argu- 
ments which quantify the creation of voids through shear: 
smaller particles preferentially fall into these voids in a 
percolation-like fashion. Therefore, it is expected that 
the larger the difference in particle sizes, the quicker this 
process will happen. While percolation rates have previ- 
ously been measured in a quasi-two-dimensional experi- 
ment [TJ, there is to date no fundamental understanding 
of the size-dependence and pressure-dependence for true 
three-dimensional flows, nor is it known in which regimes 
kinetic sieving is the dominant effect. 

We investigate the mixing and subsequent re- 
segregation of a granular material initially configured so 
that a layer of small particles is placed above an equal 
volume of large particles within an annular shear cell. 
Under shear from the bottom plate, the small particles 
migrate to the bottom and the large particles correspond- 
ingly migrate to the top, as is expected for particles of 
otherwise identical material 2 . We measure the mixing 
and segregation rates as a function of particle size ratio 
and confining pressure and find that the mixing rate is 
consistent with kinetic-sieving models for approximately- 
hydrostatic confining pressure. However, the segregation 
rates are observed to be non-monotonic in particle size ra- 
tio, in contrast with kinetic-sieving theory, and strongly 
depend on the confining pressure. Below, we quantify 
these rates and interpret them in light of the heteroge- 
neous force-transmission properties of granular materials. 

Our experimental cell is an annulus which confines the 
particles between a top plate free to move vertically and 




FIG. 1: (a) Schematic of experimental apparatus (not to 
scale) showing initial configuration of particles within the an- 
nulus. Sample images taken at window for ds — 4 mm (dark 
particles) and oIl — 6 mm (light particles): (b) initial con- 
figuration, (c) mixed state and (d) final re-segregated state. 



a rotating bottom plate, as shown in Fig. [Ha). The ra- 
dius of the inner wall of the annulus is R = 292 mm, 
and the width of the cell is 25.4 mm. Both the top and 
bottom plates are lined with rubber to increase their fric- 
tion coefficient. The bottom plate has a rotation period 
of 20.4 seconds (frequency / = 49.0 mHz), creating a 
shear band which extends a few particle diameters [5] 
into the cell. We adjust the confining pressure of the cell 
via two techniques: weighting the top plate to increase 
the pressure, or partially suspending the top plate from 
springs to reduce the pressure on the granular aggregate. 

Each experimental run begins in an initial state con- 
sisting of a layer of small particles (mass 2 kg) over a 
layer containing an equal mass of large particles. This 
initial configuration is shown in Fig. [IJb). The small 
particles are a single size for each run, with diameter ds 
ranging from 1.5 to 5.0 mm; large particles are fixed at 
diameter dr, — 6 mm for all runs. Once shear begins, the 
small particles filter downwards through the large par- 
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FIG. 2: Sample cell height data, for r = 2/6 and P = 0.36. 
An averaging window of 0.5 sec was used to smooth the raw 
signal, (a) Cell height H(t), with values H m in, Hf, and to 
marked, (b) Magnified portion of H(t) — Hmi n showing fit 
to determine mixing timescale r m . (c) Magnified portion of 
Hi — H(t) showing fit to determine segregation timescale r s . 



tides, resulting in a mixed state such as the one shown 
in (c). Eventually, nearly all of the large particles have 
reached the top of the cell (d). This experimental pro- 
tocol allows us to examine both the mixing of small and 
large particles, and the subsequent re-segregation of the 
mixture. 

The view at the outside wall is not necessarily repre- 
sentative of the bulk behavior, particularly for mixtures 
of very different sized particles. Therefore, to measure 
the average behavior of the whole system, we monitor 
the height H(t) of the top plate, as shown in Fig. |2ja). 
When shear begins, the system initially expands due to 
Reynolds dilatancy. As small particles fill the gaps be- 
tween the large particles during mixing, the overall cell 
height quickly decreases. When the particles are well- 
mixed, the aggregate takes up the least total space and 
falls to a height H min . As re-segregation occurs, the 
system re-dilates to a final height Hf. We measure the 
timescale T m for this mixing process by fitting the func- 
tion H(t) — H m i n oc e~*/ Tm to the decrease in the height of 
the cell. As the particles begin to re-segregate, the large 
particles rise through the mix and ultimately end up in a 
layer above the small particles. During this process, we 
define a segregation timescale r s by fitting a function of 
the form H(t) — Hf cx e~( t_t °)/ rs where to is chosen to 
be after the minimum i? m i n . Representative fits for T m 
and t s are shown in Fig. |2jb,c). 

As a consequence of fixing the mass (volume) of the 
particles, the total number of particles varies with par- 



ticle size ratio r = ds/d^. Therefore, we scale the cell 
height by an appropriate mean particle diameter S such 
that H/S represents the height of the cell measured in 
particle diameters. For each (rig,^) pair, we define 5 
via the relationship 
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Using 5, we can properly compare non-dimensionalized 
mixing (f2 m ) and segregation (f2 s ) rates among runs with 
different r: 



H 

*'m,s T~7 

0JT m 



(2) 



where / is the rotation frequency of the bottom plate. 

By varying the confining pressure, we explore a regime 
both above and below an approximately "hydrostatic" 
pressure due to the weight of the particles. We scale 
the effective weight of the top plate by the weight of the 
particles, and report dimensionless P defined as: 



P : 



i p g + Mg — k Ax 



(3) 



where m p = 15.42 kg is the mass of the top plate, m g = 4 
kg is the total mass of the particles, M is the added 
compressive mass (if present), and kAx is the average 
upward force from the supporting springs (if present). 
We explore values of P from 0.25 to 1.48; typical varia- 
tion within a single run is ±8 x due to the contrac- 
tion/extension of the supporting springs. We add mass 
M = to 4.5 kg to increase compression; the smallest P 
is achieved by adjusting the length of the spring supports. 

We measure the mixing and segregation rates for six 
different particle size ratios with P = 0.36 (at least 5 
runs each) and at six different pressures for r = 2/6 and 
r = 5/6 (at least 3 runs each) . Fig. [3] depicts the mix- 
ing and segregation timescales and rates as a function of 
particle size ratio with pressure held constant. We ob- 
serve that the mixing rate f2 m decreases as particles be- 
come more similar in size (r — > 1). This corresponds to 
the expected kinetic-sieving behavior [U [6j [7] whereby 
small particles filter down through a fluctuating "sieve" 
of large particles. The smaller r is, the more likely the 
small particles are to find voids to fall into. 

In contrast, we observe that the re-segregation process 
takes longer for both small and large r, as shown in Fig. [3] 
whether measured as elapsed time or a rate scaled by S. 
A maximum segregation rate is achieved near r = 3/6: 
further reductions in the smaller particle size slow the 
rate at which the system re-segregates. 

In order to better understand this behavior, it is worth- 
while to examine how £l s depends on the confining pres- 
sure on the system at both large and small r. As shown 
in Fig. [4J increasing P decreases f2 m , and the pressure 
affects contrasting particle sizes (low r) more strongly 
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FIG. 3: (a) Mixing timescale r m (o, left axis) and segregation 
timescale r s (A, right axis), as function of particle size ratio r. 
(b) Mixing rate fl m (o, left axis) and segregation rate Q s (A, 
right axis), as function of particle size ratio r. The error bars 
represent the standard error among at least 5 independent 
measurements. All data is collected at P = 0.36. 



FIG. 4: (a) Mixing rate Q m as a function of scaled pressure 
for r — 2/6 (V) and r = 5/6 (□). (b) Segregation rate Q, s as 
a function of scaled pressure for r — 2/6 (T) and r — 5/6 (■). 
The error bars represent the standard error among at least 3 
independent measurements. Dashed line is P = 0.36, which 
coincides with data from Fig. [3] 



than similar particle sizes. For fi s at low r, this effect 
is even more pronounced: a five- fold increase in pressure 
decreases the segregation rate by a factor of 100. This 
strong suppression of segregation with pressure causes an 
inversion in the r-dependence for P > 0.5. Pressure has 
little effect on either rate as r — > 1 . 

These results provide three effects in need of explana- 
tion: (1) segregation rates display non-monotonic depen- 
dence on particle size ratio, (2) contrasting particle sizes 
are much more sensitive to pressure than similar particle 
sizes, and (3) mixing rates are much faster than segrega- 
tion rates. The decrease in the segregation rate for small 
r is particularly notable, since it is inconsistent with the 
predictions of kinetic sieving. The pressure-sensitivity of 
the system suggests looking at force chains [9j [10] as an 
important factor for all three effects. 

In simulations of granular materials in two dimen- 
sions, it is observed that force chains preferentially form 
through the larger particles as size ratio r decreases 
[Til H ECU H3; this is likely related to the large par- 
ticles' enhanced number of contacts. This unequal parti- 
tioning of force chains between large and small particles, 



were it to also be present for three-dimensional granu- 
lar materials, could account for the first two effects. For 
small r, the presence of a large-particle-dominated force 
chain network at larger pressures could make it difficult 
for small particles to rearrange, thus slowing the segre- 
gation rate. As the particles become more similar in size 
(increasing r), such an imbalance would be smaller in 
magnitude. 

Another factor that could influence the anomalously 
low segregation rate of the system for small r is the ob- 
served increase in packing fraction for mixtures of dis- 
similar particle sizes |15] . If the experiments at low r are 
denser, then they have less void space and this could slow 
their re-segregation. 

We also observe a lack of reciprocity in the mixing 
and segregation mechanisms: a small particle falling 
though a mixture of mostly large particles (mixing) does 
not progress at the same rate as a large particle rising 
through a mixture of both large and small particles (seg- 
regation). Not only is there an approximately 10 x dif- 
ference in the associated rates (see Fig.|3|, but the segre- 
gation rates are much more pressure-dependent than the 
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mixing rates. The mixing process can be more clearly 
associated with the void-filling mechanisms of kinetic 
sieving, which are apparently not strongly influenced by 
pressure. However, the segregation process requires large 
particles rising (called "squeeze expulsion" by |5 ) , which 
cannot be described by void filling. 

These experiments highlight the fact that granular seg- 
regation provides a sensitive probe of how both the void 
space and the stress transmission influences the dynam- 
ics of the system. The pressure-dependency of the re- 
sults suggest that volume-based descriptors of the state 
of granular systems |16] should be supplemented by in- 
formation on the stresses [T71 [T8J [HJ HH [H]. In the 
experiments described here, we are unable to measure ei- 
ther the void distributions or the force distributions for 
large- P and small- P cases, so we cannot disentangle the 
two effects. While local free volume distributions have re- 
cently been measured in three-dimensional systems [22] , 
little is yet known how such distributions are affected by 
pressure or shear. An improved understanding of the in- 
terplay between pressure and volume state variables will 
improve models of segregation. 
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